Birkhoff's theorem fails to save MOND from non-local physics 
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We investigate the consequences of Birkhoff's theorem in general relativity (GR) and in Modified 
Newtonian dynamics (MOND). We study, in particular, the system of a finite-mass test particle 
inside a spherical shell. In both GR and MOND, we find non- vanishing acceleration for that test 
particle. The direction of the acceleration is such that it pushes the test particle toward the center 
of the shell. In GR, the acceleration is found analytically in the case of a small test mass with a 
small displacement from the center of the shell. In MOND, the acceleration is found analytically in 
the limit of large test mass and small displacement, and a comparison to numerical values is made. 
Numerical simulations are done for more general cases with parameters that mimic the system of 
a galaxy in a cluster. In GR, the acceleration is highly suppressed, and physically insignificant. In 
MOND, on the contrary, the acceleration of the point particle can be a significant fraction of the 
field just outside of the spherical shell. 
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INTRODUCTION 



Birkhoff's theorem (BT) states that in General Rel- 
ativity (GR) any spherically symmetric solution of the 
Einstein held equations in vacuum must be stationary 
and asymptotically flat. As a consequence the metric ex- 
terior to a spherically symmetric mass distribution must 
be a Schwarzschild metric. One important corollary of 
the theorem is that the metric inside a spherical shell 
(or inside the innermost of a sequence of concentric such 
sheels) is the Minkowski metric. 

Birkhoff's Theorem is the generalization from Newto- 
nian theory to GR of Gauss' Law for gravity. As a con- 
sequence of Gauss' Law, outside a spherically symmet- 
ric mass distribution, the Newtonian gravitational field 
is that of a point mass at the center of the distribu- 
tion. Meanwhile, the gravitational field anywhere inside 
any spherical mass shell vanishes, In the context of the 
inverse-square- law of Newtonian gravity, these results are 
easily understood. For example, at a point inside a spher- 
ical shell, the force on a test particle from any thin ring 
on the shell, is precisely balanced by the force due to a 
ring subtending the same angle but directly opposite the 
first ring. Even if the test particle is not at the center. 
This is because the decrease in the gravitational force due 
to one ring being farther away from the test particle, is 
precisely balanced by the increase in area (and so mass) 
of that farther-away ring. 

One can think of Birkhoff's theorem as describing the 
motion of a zero-mass test particle in the presence of 
a spherically symmetric mass distribution. A zero-mass 
test particle does not spoil the spherical symmetry. How- 
ever, in realistic situations, the "test particles" probing a 



gravitational field have a finite mass. Often, this mass is 
not small at all, as in the case of a large galaxy inside a 
cluster. Therefore, unless the test particle is at the center 
of the distribution, its presence perturbs the system from 
spherical symmetry, spoiling the assumptions underlying 
BT. This paper addresses the impact this violation of 
spherical symmetry has both in GR and in theories that 
implement the ideas of Modified Newtonian Dynamics 
(MOND). 

MOND [3, 0, 0] , the reader will recall, is an alternative 
to dark matter in which the missing gravity problem of 
galaixes is solved by altering the gravitational force law, 
rather than by introducing new unseen forms of matter. 

In both GR and MOND, there are reasons to suspect 
that the force on a test particle inside a spherical mass 
shell may not actually vanish. General relativity is a met- 
ric theory. Massive bodies change the geometry around 
them. A massive test particle distorts the space around it 
and thus a " spherical shell" ceases to be a spherical shell 
unless the test particle is at its center. This distortion 
effect from test particle will appear in any modification 
of GR that preserves its metric nature - such as DGP, 
TeVeS d, Generalized Einstein Aether (GEA) and f(R). 

In the case of MOND, which is not itself a metric 
theory (although its covariant implementations, such as 
TeVeS and GEA are), the inverse-square law which is 
crucial for the vanishing acceleration of the test particle 
inside the shell (or at least for our intuitive understanding 
of why it vanishes) does not exist. MOND was proposed 
to explain galactic rotation curves, and therefore is char- 
acterized by a gravitational field scaling as r _1 instead of 
r~ 2 at large distances from a bounded mass distribution. 

In what follows, we consider the mass distribution 
shown in Fig. Q] - a spherical shell of mass M s plus an 
interior off-center point mass, Mk, of much lower mass. 
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FIG. 1: Mass distribution: We consider a point particle inside 
a large spherical shell. R denotes the radius of the shell, and 
rfc denotes the amount of displacement of the particle from 
the center of the shell. 



We may consider this distribution to be a very simpli- 
fied version of for example a galaxy in a cluster, or a 
clump of matter - star, globular cluster, clump of dark 
matter - inside a galaxy. For defmiteness wc will the to- 
tal mass of the system to be typical of a galaxy cluster, 
M total = 2.0 x 1O 14 M , and take the radius of the shell 
to be fixed at a typical cluster raius, R = 2.0Mpc. We 
now proceed to investigate the acceleration of the point 
mass in both General Relativity and MOND. 
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Here f a — x a — x and r ao = x a — X\,. x is the spatial 
variable in the PPN coordinate frame, and x a is location 
of object a. This equation describes the free-fall of each 
point in the system. In order to describe the accelera- 
tion of the test particle in our system, we should remove 
from equation ((2]) the terms that are associated with the 
acceleration of the shell (which is much smaller than the 
acceleration of the point particle). The acceleration is 
then 
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In obtaining the final result, we have replaced the sum 
by an integral over the shell which we were able to per- 
form. Notice that there is a singularity as r k — > R. This 
singularity is not real. It appears because PPN is not 
valid as r k is very close to the shell. 

Defining d = r k /R, R Sc hs = 2GM S and g Ns - ^ 
equation ^ can be written as 
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GENERAL RELATIVITY 

In Newtonian mechanics, a spherical shell is a well- 
defined object. However, in GR, since particles change 
the geometry around them, the simple definition of a 
sphere as a set of point equidistant from some common 
point (the center) is lost in the presence of a perturbing 
mass. We define a spherical shell using Parametrized 
Post-Newtonian (PPN) coordinates. We identify a point, 
the center, such that the distance from the center in PPN 
coordinates is the same for each point on the shell. We 
call this distance R the radius of the shell. 

We can calculate the acceleration of the test particle 
in GR from the Einstein-Infeld-Hoffman (EIH) equation 
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z, (3) 



where the last line is an expansion for small d. 

The gravitational acceleration clearly does not vanish 
- it points toward the center. A particle inside a spherical 
shell is therefore attracted toward the center to restore 
spherical symmetry, where the acceleration happily does 
vanish. However, since M k <C M s , the acceleration is 
very small compared to, say, the Newtonian acceleration 
just outside the shell, gNs- It is suppressed by both the 
ratio of the test mass to the shell mass, and, more impor- 
tantly, by the ration of the Schwarzschild radius Rschs of 
the shell to its physical radius. This suppression factor 
is huge for typcial astrophysical systems, 9.7 x 10~ 6 for 
our choice of mass and radius. This suppression factor 
makes the acceleration physically insignificant in GR. 
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PERTURB ATIVE CALCULATION 



In 1984 Bekenstein and Milgrom introduced the La- 
grangian formulation of Modified Newtonian dynamics 
(MOND) 0. The field equation of MOND is derived 
from the Lagrangian 



L = — / d 3 r \ + (8ttG) 



(4) 



where ip is the gravitational potential. T{y 2 ), with y = 
is an arbitrary universal function, that together 
with ao, the characteristic scale of MOND, specifies the 
theory. Varing L with respect to tp yields a modified 
Possion equation: 



V ■ \jj,(\Vip\/a )Vip] = 4irGp(x) 
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where fx(y) = J-'(y 2 ). fx(y) must approach 1 as \y\ ^> 1 
and | y\ as \y\ <C 1, in order that the field scale as ^ 
near a spherical mass distribution (the usual Newtonian 
result) and as - far from the mass distribution to explain 
flat galaxy rotation curves. 

One the widely used form of p is 
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(However see Q for different form of ji function.) The 
value of ao is then given by phcnomcnological fit. We will 
adapt the value derived by Begeman et al [8| in the study 
of external galaxies with high quality rotation curves. 



a = 1.2 x l(T lu m/s 
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Again we consider the mass distribution in Fig. [T] The 
center of the spherical shell is chosen to be the center 
of coordinates. Because of the axial symmetry of the 
configuration, we put the point particle on the z-axis. 

We assume that the radius of the shell is large enough 
that the gravitational acceleration goes deep into the 
MOND regime. For a bounded mass distribution of total 
mass M, we define a transition radius 



Rt = y/GM/ 
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Rt indicates a point at which the Newtonian field approx- 
imately equals ao, and this is about the point at which 
the field switches from Newtonian 1/r 2 to MOND's 1/r. 
Thus we take R Rt- For the physical parameters we 
are considering, R/R t = 4.2. 

Non-linearity in the modified Possion equation negates 
the superposition principle, and it makes an analytic cal- 
culation of the field a non-trivial task. We present the 
perturbative solution for the case M k 3> M s , and in more 
general cases, we resort to numerical calculations. 



We first consider the mass of the particle inside the 
shell to be much larger than the mass of the shell. In this 
case, the density distribution of the spherical shell can 
be treated as an aspherical perturbation on a spherical 
system - the point mass - and a perturbative solution to 
the modified Possion equation is possible. Let the density 
distribution of the point particle and the spherical shell 
be pk and p s heii , respectively. Then the unperturbed field 
■00 satisfies the modified Poisson equation: 

V-[/i(|V^ |/ao)VV'o]=47rGp fc . (9) 

The exact solution for Vf/'o can be found by applying 
Gauss's theorem. In terms of the quantity u = aQ f^£ p ; 
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The expansion in small u is valid for \f — r%\ 3> y G £* k ■ 
Knowing VV>0j one can find the first-order perturbation 
equation for ipi using 



V- [|V^o|/a ((V^i -eb)eb + V^i)] = 4irGp aheU 



(11) 



where e"o is a unit vector pointing in the direction of 
Vt/V ln the deep MOND regime where /i(|V^o|/ a o) ~ 
\\7ipo\/ao, and in the limit R 3> Tk, one can expand ipi 
just inside and outside of the shell in powers of The 
first order aspherical contribution to the potential is 

<■' = <; . . (i 2 ) 
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The net force on the shell points toward positive z, having 
the effect of restoring the spherical symmetry by moving 
the shell so that the it is centered on the point mass. The 
first two terms result from tpo and the last term comes 
from ipi . We note that the first term, which is dominant 
over other terms, is independent of the specific choice 
of p function. It is a consequence solely of the MOND 
condition: p{y) — > y for y <C 1. 
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The Lagrangian ([4} from which the field equation ([5]) 
is derived is invariant under spacetime translations and 
spatial rotations. Energy-momentum and angular mo- 
mentum conservation therefore hold for an isolated sys- 
tem in MOND. In particular, Newton's third law of ac- 
tion and reaction is valid. Hence, the acceleration ak of 
the point particle in the frame of the shell can be found 
from equation (fl~3|) . Writing Rk = y/GMk/ao, d = rk/R, 
and expanding in a Taylor series, 

a<3 



+ 



As expected, at approaches zero as d — + 0. For d ^ 0, 
ak points toward the center of the spherical shell. As 
the point particle is displaced, the system reacts to re- 
store the spherical symmetry. We see here that when 
the system as a whole is not spherically symmetric, the 
acceleration on a test particle in a spherical shell is non- 
vanishing. 

The suppression of the acceleration compared to char- 
acteristic accelerations in the system is much milder in 
MOND than that in GR. In MOND, the scale of the 
acceleration is do, but suppressed by Rt/R, whereas in 
general relativity, it the scale is gjv, but suppressed by 
Rsch I R- Since Rt S> Rsch in general and since clq 3> gw 
for a system deep in the MOND regime, the scale acceler- 
ation is much larger in MOND than in GR. For M to tai = 
2.0 x 1O 14 M and Rsch — 2GM to tai, the ratio of accelera- 
tions is approximately (ao x Rt)/(gN x Rsch) = 4.3 x 10 5 . 

NUMERICAL CALCULATION 

To find the acceleration of the point particle in a more 
general case, we adapt the numerical scheme developed 
by Milgrom [10| . In this scheme, the field U is defined by 



U = /z(|V^|/a )W> 



(15) 



So long as fi is monotonic, one can invert the relation 
in equation (|15p and write V?/> in terms of U. For the 
specific choice of fj, given by equation ([5]), 
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The field U, then, satisfies the set of differential equations 

V • U = AirGp 
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FIG. 2: This figure plots the acceleration on the point par- 
ticle as a function of d — rk/R- They all points toward the 
center of the spherical shell. The numerical solutions and the 
corresponding perturbation solution is shown. 



The discretization and numerical calculation on the lat- 
tice is done in terms of U. Because of the one-to-one cor- 
respondence between and U, one can find from 
U. 

An initial ansatz for U in the numerical solution is 
given by solving 



V ■ Ui = AirGp 



(18) 



The field Ui has the correct divergence, but not the 
correct curl. The code iterates to make the curl van- 
ish at each vertex of the lattice. Ui also serves as the 
boundary condition for the numerical solution. Since 
V • 0- U~l) = 0, U and U differ only by the curl field. It 
can be shown [9( that the curl field for a bounded mass 
distribution vanishes at least as fast as ~ -4r- Then, as- 
suming that the physical size of the lattice is large com- 
pared to the mass distribution, U — ► Ui on the boundary. 

We implement Milgrom's algorithm on a spherical lat- 
tice. The point particle is placed at the center of the 
lattice, and the center of the spherical shell is displaced 
from the center of the lattice by — r^. The total number 
of angular grid lines is denoted by L, and the total num- 
ber radial grids is fixed at L/5. To meet the boundary 
condition, the radius of the outermost shell is set to be 
10(LR. 

In the case Mk 3> M s , a comparison between the nu- 
merical results and the perturbative solution is possible. 
We make the comparison between the solutions with fixed 
values of M s /M k = 0.01 and Rk/R= 0.25. Fig. [2]shows 
that these agree to better than one percent. On the 
plot shown, the percentage difference varies from 0.1% 
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FIG. 3: This figure shows the dependence of the numerical 
solution on the number of lattice sites. The value of d is fixed 
at d — 0.02. The dotted line indicates the value from the 
perturbation solution. 



to 0.9%. The discrepancy between the values increases 
as d = rk/R approaches 1, where the perturbative expan- 
sion becomes less reliable. Fig. [3] shows how numerical 
values depend on the number of lattice points. As ex- 
pected, the values approach the perturbative solution as 
L increases. In the following study, we use the lattice size 
of L = 300 and L = 600. 

Fig. 2]and Fig- E]represent the results of the numerical 
calculations. The mass dependence of the acceleration of 
the point particle is plotted in Fig. 0] and the position de- 
pendence is plotted in Fig^O We see clearly that there is 
non-vanishing acceleration directed toward the center of 
the spherical shell. In Fig. [5J we see that the acceleration 
vanishes in both the Mk — ► limit and as M s — ► 0. This 
is expected since these are two limits in which spherical 
symmetry is recovered. For the three curves plotted in 
Fig. 31 the peak value occur around Mk/Mtotal ~ 0.15. 
The position of the peak has a slight position dependence 
- as the particle gets closer to the shell (large d) , the peak 
occurs at smaller Mk/Mtotal ratio. 

We note one curious feature from Fig. 0] The astro- 
nomically interesting region in Fig. |4]is where Mk/M 8 <C 
1. In this regime, the acceleration of the point particle 
is a very sharp function of its mass. This implies that 
galaxies with slight mass differences might experience 
quite different acceleration in theoriest that implement 
the MOND limit. 

The distance dependence of the acceleration is shown 
in Fig. El Again, the acceleration of the particle van- 
ishes when the particle is near the center (d — > 0), 
and increases monotonically outward. Near the shell 
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FIG. 4: This figure shows the acceleration of the point par- 
ticle as a function of Mk/Mtotal for particles at a variety of 
positions. The acceleration points toward the center of the 
shell. The numerical results are for L=300. 



(d ~ 1), the value for Mk/Mtotal — 0.01 exceed that for 
Mk/Mtotai — 0.15. This is because the peak mentioned 
above occurs at lower value of Mk/Mtotal for larger value 
of d. 

Unlike for GR, in MOND the acceleration of the point 
particle is a significant fraction of <?at. Especially when 
the particle is close to the shell, the acceleration can be 
larger than . 



CONCLUDING REMARKS AND IMPLICATIONS 

In both general relativity and MOND, the accelera- 
tion of a massive test particle inside a spherical mass 
shell vanishes only when the particle is at the center of 
the shell. When the particle is displaced from the center, 
the particle experiences a force toward the center of the 
shell, the direction in which spherical symmetry would 
be restored. The magnitude of the acceleraion in GR is 
not physically significant in most (or probably all) astro- 
physical situations of interest, since it is suppressed by 
the ration of the Schwarzschild radius of the shell to the 
size of the shell. 

In MOND, on the contrary, the acceleration of the 
point particle is a significant fraction of the surface grav- 
ity just outside the shell. This is despite the fact that 
MOND has a Birkhoff or Gauss-like theorem which im- 
plies that the potential inside an isolated spherical mass 
shell is constant. However, this theorem operates in 



6 




FIG. 5: This figure shows the acceleration on the point par- 
ticle as a function of position for fixed particle-mass to total- 
mass ratio. Results are plotted for two different mass ratio. 
The acceleration points toward the center of the shell. Simu- 
lations are performed on the lattice with L=600. 



the absence of the usual explanation for Gauss' Law in 
Newtonian gravity or classical electrostatics - the bal- 
ance between the r 2 growth of surface areas and the r~ 2 
force law. Thus, when the very particular conditions of 
Birkhoff's theorem are broken even a little, the gravi- 
ational force re-emerges at considerable strength. The 
characteristic MOND field for the values of parameters 

VGMtataK 



From Fig. [4] and Fig. [5j we can see that the acceleration 
inside the spherical shell can be significant at this scale. 

In MOND, the applicability of Birkhoff's theorem is 
very limited, and in computations one needs to consider 
not only the local mass distribution but also the back- 
ground mass distribution. 
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